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Z: ■ ABSTRACT 

Oh, Different approaches to the infrared emission spectrum of molecules are dis- 

O . cussed with a view to understanding the spectral profile of the Unidentified 

+^ ' 
C/3 ' Infrared Bands. Normal mode analysis, which delivers the integrated band 

intensities, is found to be inadequate: in particular, it does not account for 

the band widths. A more accurate procedure requires the computation of the 

\^ i energy contents of the vibration modes as well as the variations of the total 



electric dipole moment of the excited molecule. These quantities reflect the 
anharmonicity of atomic motions in molecules above the ground state, as a 



f~^ , consequence of which neighboring modes interact and exchange energy. This 

has far reaching effects: a) the initial distribution of energy among modes, 
just after excitation, is far from uniform; b) energy is redistributed during ra- 

^^ ' diative relaxation, so the energy ultimately radiated by a given mode differs 

j^ , from the energy that was initially imparted to it; c) vibration lines are broad- 

ened and acquire wings; d) as the molecule grows in size, line wings overlap 
and form plateaus underlying regions of high line concentration; e) this is ac- 
companied by an increase of the number of phonons, vibrations involving the 
whole molecular structure, whose accumulation ultimately contributes to the 
far infrared continuum; f) the emission spectrum differs considerably from the 
normal mode spectrum and depends on the excitation process. 

This is illustrated by computations on molecules of increasing size: 6, 
101 and 493 atoms. A tentative UIB-like synthetic spectrum is obtained by 
summing up the spectra of 21 CHONS structures, computed in the same 
way. Also discussed is the relative contributions of the molecule to band and 
continuum emission, as a function of its size. 

Key words: astrochemistry — ISM:lines and bands — dust. 
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1 INTRODUCTION 

The quest for a realistic model of the structure, composition and excitation 
of the carriers of UIBs (Unidentified Infrared Bands) has been going on for 
more than three decades now. An essential requirement is, for such a model, to 
correctly reproduce the spectral features of UIB emission, which, by now, are 
very well documented in a large variety of astronomical sites. Despite some 
variations from object to object, there are enough similarities that it has been 
possible to partition most of the observed spectra into a small number of 
classes (see, for instance, Peeters et al. 2002, 2004; van Diedenhoven et al. 
2004), thus allowing systematic comparisons to be made with laboratory or 
theoretical spectra. 

One universal characteristic of UIBs is that they are relatively wide, and 
continuous even at a resolution of 1500; moreover, "plateaus" appear beneath 
the stronger bands, between 6 and 8 jiva and between 11 and 13 /im (see, for 
instance, Hony et al. 2001, Kwok and Zhang 2011, etc.); also, an underlying 
continuum may be seen to increase steadily from 6 /im upwards (see Smith et 
al. 2007). The present paper focuses on these particular aspects of UIBs. 
Most present day candidate carriers are big molecules of widely different sizes 
(see Tielens, 2008), the size being limited, of course, by experimental or com- 
putational constraints. While the experimental production and spectral anal- 
ysis of single macromolecules have proven to be daunting, modern computing 
capabilities make it easier to use theory or chemical simulation. The most 
usual computational way of determining the IR (InfraRed) spectrum of a 
given structure is to perform a Normal Mode Analysis (see Wilson et al. 
1955). Once a composition and a structure have been chosen for the model 
carrier, the first step of this procedure is to compute the ground state in vac- 
uum at K, i.e. the state of lowest potential energy. A classical mathematical 
transformation then delivers all the different modes of atomic vibration of in- 
finitesimal amplitude, characteristic of the grain. In the usual approximation, 
such a vibration corresponds to the motion of a massive point about the apex 
of a parabolic potential well. For a grain made up of N chemically bound 
atoms, there are 3N-6 such modes, each characterized by its frequency and 
the directions and velocities of the excursion of each atom from its equilib- 
rium position (the apex), along the 3 space coordinates. From these, and the 
atomic charge distribution, the dipole moment can be deduced for each mode. 
Further treatment yields the mode IR intensity (or line intensity, or integrated 
hand intensity), which is proportional to the absorbance at the corresponding 
frequency, the quantity commonly measured in the laboratory. 
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While this procedure is higly instructive, it suffers, for our present purposes, 
from at least three defects: 

a) The spectral mode intensity by itself is not sufficient to predict emission 
intensity at the corresponding frequency. One also needs to know the extent to 
which this mode is kinetically excited; this, in turn, depends on the excitation 
process and on the coupling between modes. 

b) As a result of excitation, emission lines may be shifted from absorption 
frequencies. 

c) By definition, the normal mode line intensity is a (5 function and is therefore 
devoid of any width. By contrast, emission lines are broadened by anharmonic- 
ity of intramolecular atomic motions. 

As early as 1992, Brenner and Barker (|1992p warned that "unambiguous iden- 
tification of the IS (Interstellar) emitter cannot be established just on the 
basis of matching the emission frequencies to laboratory absorption spectra" . 
Brenner and Barker themselves computed the emission of excited benzene 
and naphtalene based on a) the ergodic assumption (statistical distribution 
of energy), b) Einstein's spontaneous emission coefficients, as deduced from 
measured (relative) absorption coefficients at a few characteristic spectral fre- 
quencies of the two molecules (see also AUamandola et al. 1989). Intrinsic 
bandwidth and underlying emission continuum were assumed. 
The emission bandwidth may be the result of a number of processes acting 
together (e.g. Cook and Saykally 1998, Papoular 1999, Pech et al. 2002). Some 
width components depend on the energy content of the mode, and, therefore, 
on the excitation energy and the grain size. An essential broadening cause is 
the continuous exchange of energy between modes coupled by anharmonicity. 
In this respect, comparison with laboratory measurements may be misleading. 
For, in the laboratory, the detection and spectroscopy of molecular bands usu- 
ally require much higher densities and temperatures of molecules than found 
in interstellar space (e.g. Nemes et al. 1994, Frum et al. 1991), in which case 
pressure broadening may far exceed other broadening effects (see, e.g., Kuhn 
1971, Peach 1981). Another common broadening effect in the laboratory is 
molecular rotation (e.g. Pirali et al. 2009), which must be quite weak for large 
molecules in thin and cold astronomical environments. As for anharmonicity, 
it is not easy to determine, as it varies from vibration mode to vibration mode 
and from molecule to molecule (see Brenner and Barker 1992). For want of 
better knowledge of the broadening mechanisms in action, an arbitrary broad- 
ening is often convolved with computed normal mode intensities (e.g. Boersma 
a2010, b2011). An ad hoc frequency shift may also be necessary to fit models 
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with observations (e.g. Verstraete et al. 2001). The present state of the art 
in this domain is clearly summarized by Bauschlicher et al. (|2010p in their 
Appendix. 

Here, these issues are addressed by using computational chemistry codes to 
deduce the emission intensity from, the variations of the dipole moment of the 
particle after it has been excited. In this way, the most important broaden- 
ing processes are included, as well as the particulars of the excitation and 
relaxation processes, and the particle size. This method may be applied to 
vibrational emission following any type of excitation mechanism. No a priori 
assumption (such as ergodicity or thermal equilibrium) is made regarding the 
relaxation process. The present work seeks to show why and how the dis- 
crete vibrational spectra of molecules evolve into continuous band spectra, 
as observed. It takes advantage of the huge progress that occured in digi- 
tal computation chemistry during the last decades, which makes it possible, 
today, to apply long known theoretical notions to relatively large molecules. 
In a way, this work is also a complement, or a twin, to a previous paper 
describing in detail the excitation by H atom impact, and the ensuing rate 
of vibrational emission, regardless of the spectral profile (chemiluminescence; 
Papoular 2011a). 

To begin with. Sec. 2 describes three approaches to the spectrum, in order of 
increasing accuracy, as applied to a very simple, planar, molecule, ethylene 
C2II4. This highlights the strengths and deficiencies of each approach. Atten- 
tion is drawn to the finding that internal thermal equilibrium (or ergodicity) 
does not usually apply after excitation of the molecules, contrary to a common 
simplifying assumption (see, for instance, Brenner and Barker 1992, Mulas et 
al. 2006). 

The third, and more accurate, procedure, based on time analysis of the dipole 
moment variations, is then applied, in Sec. 3, to larger molecules (101 and 
493 atoms), in which lines are broadened enough to form bands in crowded 
spectral regions. These molecules contain a variety of atoms of the "CHONS" 
family (C,II,0,N, and S), some of the most abundant in space. Section 4 
carries the endeavor a step further, by applying the dipole moment procedure 
to several more molecules of similar composition and 27 to 112 atoms in size 
each, and adding their emission spectra. As a result, plateaus, and even a 
budding continuum, become apparent. In Sec. 5, the expected evolution of 
the relative contributions of bands and continuum to emission is discussed as 
the molecules grow into grains by clustering or sticking. The main outcomes 
of this research are summarized in the Conclusion. 
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2 THREE APPROACHES TO VIBRATIONAL 
SPECTRA 

The present work is based on the use of various algorithms of eomputational 
organic chemistry, as embodied in the Hyperchem software provided by Hy- 
percube, Inc., and described in detail in their publication HC50-00-03-00, and 
cursorily, for astrophysical purposes, in Papoular pOQip . Here we use the im- 
proved version Hyper 7.5. While almost all the computations used standard 
algorithms, some required writing special purpose subroutines, a capability 
also provided by the software. Semi-empirical methods of quantum chemistry 
(PM3 and AMI) were used for molecules up to about 100 atoms in size. Be- 
yond that, the less accurate, but faster Molecular Mechanics (MM-f-) method 
was preferred. Anharmonicity is built in and tailored by comparison with 
laboratory experiments on hydrocarbon molecules. 

2.1 Normal Mode Analysis 

Some of the main points of the present argument will now be made through 
chemical "experiments" on ethylene. This was selected because it is small 
and, hence, amenable to more accurate simulation; the modes are few, far 
apart and easily distinguishable in frequency and behavior; and because the 
internal dynamics of its atomic vibrations were studied in detail in Papoular 
(|2006|) . The first step, in all approaches, is the building and optimization of 
the structure, i.e. minimizing its potential energy by automatic standard pro- 
cedure. This is found to be -533.63 kcal.mol~^(in PM3; 1 eV=23 kcal.mol^^). 
A standard NMA then delivers the integrated IR band intensity, I{i'), which 
is linked to the corresponding absorbance in cm~^, at frequency v (cm^^), by 

a(^) = 10^/M^ , (1) 

where Av is the band width (cm^^) and C is the molecular density (mol.l~^). 
The computation of / is independent of Aj/, and i> is the frequency of an ideal 
vibration of infinitesimal amplitude about the equilibrium atomic geometry. 
For ethylene, which has 12 modes, the intensity spectrum is shown in Fig. 1. 
By definition the NMA and spectral intensities are best suited to the predic- 
tion of the outcome of an absorption experiment, as this is done on samples 
which are (nearly) in their ground state. In the case of ethylene, only 3 of 
the modes are notably "IR active": 2 (antisymmetric) CH stretching modes 
and 1 oop (out of plane), in phase CH vibration. Simple inspection of the 
molecular geometry shows that, in the other modes, the local dipole moments 
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Figure 1. Infrared spectral intensities of ethylene C2H4. These are proportional to the corresponding absorption coefficients. 
Note the logarithmic ordinate scale. In effect, only 3 modes have enough intensity for them to qualify as IR-active: 2 stretching 
modes, near 3.2 ^m, and the in-phasc, out-of-planc CH bending mode at 10.2 ^m. The others are "dipole-momcnt forbidden". 



resulting from the movements of different atomic groups mostly cancel each 
other. This, of course, is a precious piece of information. However, by itself, it 
is not sufficient to predict an emission spectrum. Obviously, this also depends 
on the energy content of each mode, which we take on next. 



2.2 Energy distribution among modes 

An isolated molecule in space can aquire energy by, for instance, absorbing 
a photon, or by one of its dangling (C-) bonds capturing an H atom, thus 
increasing its potential energy by an amount equal to the binding energy 
(Papoular 2011a). In the following, we assume that the average interval be- 
tween two successive excitation events is much longer than the IR lifetime 
of the molecule. So, the radiative power of a molecule is the product of the 
frequency of exciting events and the energy radiated by the molecule upon 
radiative relaxation. If the various vibration modes of the molecule were not 
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coupled, all the energy deposited somehow in a particular mode would ulti- 
mately be reradiated at the frequency of the said mode. In fact, the modes 
are coupled by anharmonicity, so the modes which radiate efficiently, i.e. at a 
high rate, are constantly "refueled" by those which radiate more slowly, until 
all the energy deposited by the exciting agent has been dissipated. Ultimately, 
therefore, the former modes will radiate more energy than the latter. Further 
analysis (see Appendix) shows that the total energy radiated by mode i/ upon 
relaxation is proportional to 

E*{iy)^u^I{iy)Eo{iy), (2) 

where Eq{v) is the initial energy carried by the said mode, before radiative 
relaxation starts; this is assuming that the energy deposited by the exciting 
event is redistributed among the various modes at a much higher rate than 
the rate of radiative emission. Indeed this intramolecular vibrational redistri- 
bution (IVR) is usually completed within a few nanoseconds, much shorter 
than the radiative lifetime (up to 1 s). 

We therefore seek, now, to determine how the energy deposited in a carrier by 
the exciting agent is initially shared among the various modes. This partition, 
E{i>), where the subscript was dropped, has no direct link with the spectral 
intensities and may depend on the excitation process. Here, the focus is on 
excitation by H atom capture by a carbon dangling bond supposed to preexist 
at the surface of a hydrocarbon grain (Papoular 2011a). In this process, a 
free H atom approaches the grain closely enough to feel the attraction of a 
carbon atom whose bonds are incompletely occupied. Under favorable initial 
conditions, it falls into the corresponding potential well. It was shown that 
this capture reorganizes the distribution of the constitutive atoms in space so 
as to lower the ground state of the newly completed structure by an amount 
equal to the CH bond energy {Ei, ~4.5 eV). This is equivalent to depositing 
this energy in this new structure, an energy which is immediately and evenly 
divided into kinetic energy and potential energy (relative to the ground state 
of the completed structure). The kinetic energy itself, initially localized in 
the CH stretching vibration, is then redistributed among this and other (not 
necessarily all) molecular vibrations. This could not occur were it not for the 
couplings between modes. While normal modes are uncoupled by definition, 
they become coupled as soon as some energy is deposited in the molecule, 
because of anharmonicity, which sets in above the ground state. 
Note that, in this process, by contrast with UV photon excitation, the molecule 
remains in the electronic ground state even though atomic vibrations are ex- 
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cited. To simulate this process on the computer, we start with an optimized 
structure (here, ethylene) at rest, with all bonds occupied. One of the 4 CH 
bonds is then rotated out of the molecular plane by several degrees and ex- 
tended from its original length by a few hundredths of A so as to increase the 
molecular potential energy by, say, E, which is tantamount to depositing this 
energy in the molecule. The molecule is then set free to perform a molecular 
dynamics run in vacuum, during which all atomic motions are monitored on 
the computer screen and recorded, as well as the total kinetic and potential 
energies, as a function of time. 

Initially, as expected, the dominant motion is the CH stretch accompanied 
by excursions out of the molecular plane. Soon, however, the other atoms 
are dragged along: inspection and comparison of the variations of, say, the 
bond lengths immediately reveals that energy constantly flows from bond 
to bond, thus increasing the amount of energy shared by other bonds. This 
process translates into a continuously decreasing standard deviation, a, of each 
and every variable, like the total kinetic energy of the system, for instance. 
After some time (the relaxation time), a dynamic equilibrium sets in, which is 
characterized by constant average relative standard deviations. The relative 

1/9 

standard deviation of the energy, as/E is roughly A^^ , where Nm is the 
number of effectively excited normal modes (see Reif 1984). The relaxation 
time is highly variable, depending on the grain structure and size, and on its 
initial perturbation. It may be as long as hundreds of picoseconds. This is still 
much shorter than the minimum radiative lifetime (^ 1 ms). 
Once the dynamic equilibrium is reached, a new run is launched, during which 
all variables of interest are monitored. Here, we focus on the instantaneous 
total kinetic energy. Its FT (Fourier Transform) consists of peaks at twice the 
frequencies of those of the intensity spectrum, i.e. those of the normal modes 
(Sec. 2.1). What may vary (moderately) from run to run is the peak heights. 
These peaks are artificially broadened by the finite length, A t, of the run: 
their width at half maximum (FWHM, cm^^) is ^ (cAt)~^. The longer the 
run, the more representative is the spectrum, of the average energy content 
of the modes: 10 ps will do. 

Figure 2 illustrates the procedure as applied to the ethylene molecule, after 
it was energized as explained above: one of the CH bonds was rotated out of 
the molecular plane by 20 deg and extended from its original length, 1.086 
A , to 1.15 A so as to increase the molecular potential energy by _E ^ 5 
kcal/mol. After allowing some relaxation time, the total kinetic energy was 
recorded over a period of 10 ps. Its excursions from its average value, ^ 2.2 
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Figure 2. The curve illustrates the spectral distribution of molecular kinetic energy of ethylene after perturbation of one CH 
bond, represented here by the amplitude of the FT of the instantaneous kinetic energy variations during a 10 ps run of molecular 
dynamics on the free molecule in vacuum (sec text). The dots mark the wavelengths of the normal modes. The average kinetic 
energy is 2.2 kcal/mol (vibrational "temperature": 122 K). This spectrum is distinctly sparser than that of the IR spectral 
intensities, Fig. 1, because the normal modes are only weakly coupled. Nonetheless, this coupling gives rise to combination 
bands (intermodulation) between modes, which arc visible upon zooming, but do not coincide with any of the dots. 



kcal/mol, were extracted. The FT of the latter was delivered by a FFT (Fast 
Fourier Transform) algorithm, and its amplitude is represented in Fig. 2. This 
is a measure of the spectrum of the kinetic energy density; the scaling factor 
somewhat depends on the filter window used by the FFT algorithm, and on 
the length of the molecular dynamics run. 

Only 2 significant peaks emerge in Fig. 2: one coinciding in wavelength with 
CH stretching normal modes, near 3.2 /im, and the o.o.p. (out of plane) CH 
bending mode near 11.3 /im (cf. Fig. 1). The peak near 6.3 /im is a sub- 
harmonic of the strong stretching mode. To see if this result depends on the 
excitation process, a pristine molecule was first drowned in a thermal bath 
at 140 K so that it acquired an internal kinetic energy about equal to that 
injected by perturbing one of the CH bond (as in the previous experiment). 
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It was then freed to vibrate in vacuum, as above. The resulting kinetic energy 
spectrum is shown in Fig. 3: ahhough the stretching is stih dominant, Fig. 3 
is visibly different than Fig. 2, indicating that the spectrum does depend on 
the details of the excitation process; but there still is no energy equipartition 
between modes. Equipartition is approached only when the spectrum is taken 
on the molecule drowned within the bath, and continuously exchanging energy 
with it, which, of course, is hardly the case when the molecule is in a thin 
medium, and excited only once in a while. This confirms that, at such low 
energy densities (^ 0.2 kcal/mode), in so small a molecule, the energy dis- 
tribution in dynamic equilibrium, in vacuum, a short while after excitation, 
mainly depends on coupling between modes; and the vibration modes are only 
weakly coupled except for very close modes such as the CH stretchings. In- 
deed the records of the variations of the CH bond lengths, show continuous 
exchanges of energy between these bonds, not the others. 
These continuous exchanges translate into side bands on the wings of the 3.3 
fi peak, which can be made visible by zooming on Fig. 2 and 3. For large 
molecules, these side bands contribute to forming underlying plateaus and 
continua in spectral regions of dense line clusters. A more detailed discussion 
of anharmonicity and mode coupling effects in this molecule can be found in 
Papoular f 20061) . For our present purposes, note, in Fig. 2, the absence of the 
0.0. p., in phase, CH bending mode at 10.2 /im, which Fig. 1 shows to have 
the strongest IR intensity. This is because the asymmetry of the particular 
excitation process applied here is incompatible with the symmetry of this 
mode (all 4 CH bonds vibrating in phase), so the latter cannot be excited. 
As a consequence, if, using Fig. 1 and 2, a plot of v'^I{v)E{v) is drawn to 
represent the spectral emission of the excited molecule (eq. 2), it turns out 
that only two CH stretching modes emerge from the lot: the others either 
are too weakly excited or have weak IR intensity. However, as the /'s are 5 
functions at fixed wavelengths, this procedure runs into difficulties as soon as 
anharmonocity becomes strong enough to notably shift and/or broaden the 
spectrum of E{v) w.r.t. that of /(i^). One is then led to the third, and more 
rigorous, procedure. 



2.3 Dipole moment analysis 

Chemical simulation experiments, such as described above, highlight the im- 
portance of anharmonicity as soon as energy is injected in a molecule. On the 
other hand, normal mode intensities are only valid in the ground state, or in 
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Figure 3. As in Fig. 2, the curve illustrates the spectral distribution of molecular kinetic energy of ethylene, represented by 
the amplitude of the FT of the instantaneous kinetic energy variations during a 10 ps run of molecular dynamics on the free 
molecule in vacuum (see text). Here, however, the excitation consisted in immersing the molecule in a thermal bath for a few 
picoseconds, until it reached thermal equilibrium with about the same internal kinetic energy as in the previous experiment. 
It was then left free to vibrate in vacuum, and monitored as above, so as to obtain the energy spectrum. The dots mark the 
wavelengths of the normal modes. The spectrum is different, but not much richer, than that of Fig. 2. 



the limit of infinitely small exeursion from ground state. As a matter of fact, 
the I's are specifically defined for absorption from the ground state to the 
first excited state. Now, the electromagnetic emission of a vibrating dipole 
is known to be proportional to the square of its electric dipole moment. A 
more realistic way of obtaining the emission spectrum in response to a given 
excitation is, therefore, to turn to the time variations of the electric dipole 
moment of our excited structure. Let this be M{t), < M{t) > its average and 
M its variable part. Each frequency component, m^, of M corresponds to one 
normal mode of frequency i>. As shown in the Appendix, it turns out that 



I{v) = 2.6vml 



(3) 



where / is in km/mol, v in cm ^ and m in Debye. It follows from eq. 2 that 
the radiated energy scales like 
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Figure 4. The bottom graph represents the spectrum of the energy or power radiated by the molecule after excitation by bond 
extension, as in Fig. 2. This was obtained by monitoring the total electric dipolc moment of the molecule and applying eq. 4. 
For comparison, the middle graph represents the normal modes of Fig.l (linear ordinates here), and the upper curve represents 
the kinetic energy spectrum of the excited molecule, as in Fig. 2. Most normal mode lines are conspicuously absent in emission 
, because the corresponding vibration modes are either IR-inactive or not excited. 



E* = v^I{v)E{v) ^?:,lQ-^v''miE{v) , 



(4) 



where E{v), in eV, is the initial energy distribution, as obtained in Sec. 2.2. 
The first expression of E* in eq. 4 is strictly valid for harmonic motions only, 
when I{v) retains its original meaning, This is not the case of the second 
expression, for the dipole moment incorporates anharmonicity and broadening 
and shifting effects, if any. From now on, the latter expression will therefore 
be used to represent the radiation emission of the excited molecule. Now, 
chemical simulation codes also deliver m{v^t): during a molecular dynamics 
run, these codes also monitor the displacement of electrical charges through 
the structure (using the Born-Oppenheimer approximation) and compute the 
3 Cartesian components of the resulting overall dipole moment. To derive the 
squares of the individual frequency components of the latter, one computes the 
sum of the autocorrelation functions of the 3 components then takes the FT of 
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Figure 5. This is a zoom on the CH stretching region of Fig. 4. As a consequence of excitation and anharmonicity, the 
separation between the two IR-active CH-stretching hnes is nearly twice that of the corresponding normal modes, and the 
shorter wavelength emission peak is dominant. The horizontal segment, of length 0.0012 fira (bottom graph), represents the 
artificial broadening due to the finite calculation run time (30 ps). The real, physical, broadening is ^0.004 fim. The relative 
shift of ^ 0.035^m (1 %) between the middle spectrum and the other two is probably to be ascribed to minor inconsistancies 
between the two distinct algorithms which compute normal modes and molecular dynamics, respectively. 



this sum. This was done with the results of a run 10 ps in length (3333 cycles, 
0.003 ps each in length; the dipole moment being computed after each cycle) 
and the bottom graph in Fig. 4 plots the product v^m'^E{v), to represent the 
radiation energy ultimately emitted (eq. 4), as a function of wavelength. For 
comparison, the middle graph represents the normal modes of Fig.l, and the 
upper curve represents the kinetic energy spectrum of the excited molecule, 
as in Fig. 2. Most normal mode lines are conspicuously absent in emission 
, because the corresponding vibration modes are either IR-inactive or not 
excited. 

Figure 5 is a zoom on the CH stretching region of Fig. 4. It reveals more 
subtle effects: as a consequence of excitation and anharmonicity, the separa- 
tion between the two emission peaks of CH-strctching modes is twice that 
between the two corresponding normal mode peaks, and the ratio of the two 
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emission peaks is reversed. Only dipole moment analysis can reveal such, and 
other, effects which go undetected by Normal Mode analysis or Kinetic Energy 
analysis alone. 



3 FROM LINES TO BANDS 

While the simple ethylene molecule is convenient for demonstrating the pro- 
cedure and value of Dipole Moment analysis, it also hints at the fact that 
small molecules (comprising less than a few tens of atoms) cannot carry the 
UIBs. This is also borne out by ample chemical simulation and observational 
evidence (see Kwok and Sandford 2008, Kwok and Zhang 2011). Obviously, 
this is because they cannot contain enough variety of sites and atomic species 
to match the variety of vibration features encountered in UIBs. 
In molecular spectra, it is possible to identify, in the near and mid-IR, groups 
of vibration lines carried by, and characteristic of, small groups of bound atoms 
called functional groups. Each of these is composed of the same chemical 
species associated in similar clusters, which may be replicated at different 
locations in the molecule. Depending on their environment, distinct functional 
groups of the same type will give rise to neighbouring but separate lines. 
One may therefore hope to build a band by accumulating functional groups 
in slightly different environments. Small molecules lack space to carry many 
distinct groups of the same type; so, even a very large number of identical 
small molecules will not do. 

The situation improves with medium-sized (size of order 100 atoms) or large 
(of order 1000 atoms in size) molecules, as the lines crowd in each near- or 
mid-IR band. As the lines get closer, anharmonicity couplings become tighter. 
Those lines which are close enough to merge into bands are obviously carried 
by the modes of particular functional groups: CH stretchings, bendings (i.p. 
and o.o.p.), etc. The band forming process goes on with increasing molecular 
size, provided the average chemical composition and structure are preserved. 
This is precisely how (disordered) solid state bands are formed. Besides, if the 
molecule is large enough, the wings of adjacent broadened lines encroach over 
each other, to form plateaus and, ultimately, an underlying continuum. 
We therefore turn now to a more realistic model carrier, viz. a medium sized 
molecule, to monitor the spectral changes as the size increases. Figure 6 rep- 
resents an organic ( so-called CHON) molecule comprising 51 carbon atoms, 
39 hydrogen atoms, 8 oxygen atoms, 1 nitrogen and 2 sulfur atoms, for a total 
of 101 atoms. It is composed of 1 pyrene (4 pericondensed (compact) benzene 




Figure 6. Aromatic and aliphatic structure with 101 atoms of the CHONS family. As usual, C and H atoms are not marked 
explicitely. 



rings) and 2 groups of 3 catacondensed (linear) rings including pentagonal 
rings. These substructures are linked together by CH2 and O bridges. All 
three procedures described in Sec. 2 were applied. 

To simulate the capture of an H atom by a dangling (unoccupied) C bond as 
in chemilumunescence (see Introduction) , the length of one of the CH bonds of 
the molecule was extended from 1.1 to 1.2 A . This increased the total poten- 
tial energy by 104 kcal/mol (4.5 eV). The molecule was then allowed to relax 
freely in vacuum for 22 ps, during which potential energy progressively de- 
generated into kinetic energy until a dynamic equilibrium was reached, where 
both energies became nearly equal on average. The dipole moment spectrum 
was then obtained from a molecular dynamics run 10 ps in length, as in Sec. 
2. 

Figure 7 shows the normal mode spectrum (center graph, 297 lines) together 
with the energy distribution (upper graph) and the molecular dipole moment 
spectrum (bottom graph). It is observed again that many normal mode lines 
do not show up in the latter: they are either "dipole forbidden" or simply not 
mechanically set in motion. Most importantly for our present argument, many 
lines merge into bands. This is illustrated in Fig. 8, by zooming in on the CH 
stretching region: the 58 CH stretching lines coalesce into essentially 2 bands, 
near 3.24 and 3.28 /im respectively, corresponding to CH bonds in different 
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Figure 7. Each of the 3 graphs corresponds to one of the approaches (Sec. 2) to the energy radiated by the lOl-atom molecule 
sketched in Fig. 6, when excited as described in the text. Each emission peak corresponds to a group of features appearing in 
both the energy and the IR intensity graphs (cf. eq. 4): the emission spectrum is therefore sparser than both the other two. 
Note the budding continuum underlying the kinetic energy spectrum. 



environments: 0-capped pentagons and aromatics, respectively. Another ex- 
ample of this line coalescence is to be found in the CH o.o.p. bending region, 
11-12.5 /im. Thus, bands were formed, whose width cannot be due to finite 
run time: the broadening due to the latter (10 ps here) is only 3.510^'^/im. 
The cause is, rather, the fast and unceasing exchange of vibrational energy 
between close modes, which become coupled by anharmonicity as soon as en- 
ergy is deposited in the molecule; this is amply documented by computational 
experiments (see for instance Papoular 2006). Here, this is made possible by 
the larger numbers of atoms, and hence of normal modes, so the latter come 
closer to each other. 

The width of the aromatic CH stretching emission band at 3.28 /im is now 
0.016 /im. This is not much narrower than Type 2 3.3-/im astronomical band 
of Tokunaga et al. (|T99lp or van Diedenhoven p004p . but still less than half 
his Type 1 (0.04 /im). A still larger molecule was therefore built by linking 
5 copies of the lOl-atom molecule of Fig. 6 through short bridges, to give 
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Figure 8. This is a zoom of Fig. 7, on a typical spectral region, with a resolution <5A ^ 0.0035 /im. The short horizontal segment 
represents this artificial broadening due to the finite calculation run time (10 ps). The much larger, physical, line broadening 
is not artificially enforced. It is due to various effects of intrinsic anharmonicity revealed when the molecule is excited above 
the ground state. Also, close lines coalesce into bands. A related consequence is the piling up of wings of close lines, forming 
plateaus beneath the peaks between 6 and Sfim, 11-12 ^ni, etc. 



a structure of 493 atoms. Using the PM3 semi-empirical method, Normal 
Mode Analysis delivered the spectrum of 1473 line intensities shown in Fig. 
9, center graph. Such big molecular sizes are not easy to process with the 
semi-empirical algorithms designed for PC's, like Hyperchem: Fig. 9 required 
720 h of calculation time. Performing Dipole Moment Analysis with the same 
chemical simulation method proved hopeless. It is, however, possible to resort 
to the less accurate, but much faster Molecular Methods, such as MM+ (see 
Papoular 2001), which is used hereafter. 

This molecule was first excited, as above, by extending one of the CH bonds 
from 1.1 to 1.68 A , which gave it an excess potential energy equivalent to 
that which it would gain by capturing an H atom on a dangling bond (4.5 eV; 
see Sec. 2.2). After freeing it to relax in vacuum, a dynamic equilibrium was 
reached with a total internal kinetic energy of ^ 56 kcal/mol, corresponding 
to a temperature of '-^ 38 K. A Dipole Moment Analysis was then performed 
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Figure 9. Each of the 3 graphs correspond to one of the approaches (Sec. 2) to the energy radiated by a "CHONS" molecule 
carrying 493 atoms, built up by linking together, through short bridges, 5 molecules, each identical with that shown in Fig. 6, 
and excited as described in the text. The same comments can be made as on Fig. 7. 



on a molecular dynamics run of 3333 cycles, 0.003 ps each, for a total of 10 

ps, ensuring negligible line broadening due to finite run time. The resulting 

spectra are shown in Fig. 9 and 10. As a result of the larger size of the molecule, 

the line broadening and coalescence signaled in Fig. 7 and 8 are only enhanced 

here. 

Another, related, consequence of the increased molecular size is the growth 

of plateaus in regions of high normal mode density: between 3.2 and 3.4 /im, 

8.5 and 10 /xm, 11 and 12 /xm, etc. This is a result of the wings of close lines 

spreading out over each other. 

By comparison with Fig. 8, note the appearance, in Fig. 10, of new lines 

beyond 3.3 /im, due to new aliphatic sites in this molecule. This illustrates 

again how the number of close, but not coincident, IR lines can be made to 

increase in any spectral region characteristic of a functional group, so as to 

approach the shape of observed bands. 

A potential problem with the Molecular Methods of chemical simulation, when 

it comes to absolute values, is the frequency error intrinsic to the code's use 
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Figure 10. This is a zoom of Fig. 9, on a typical spectral region, with a resolution SX ~ 0.0035 ^m, represented by the short 
horizontal segment. The much larger, physical, line broadening is not artificially enforced. It is due to various effects of intrinsic 
anharnionicity revealed when the molecule is excited above the ground state. Also, close lines coalesce into bands. A related 
consequence is the piling up of wings of close lines, forming plateaus beneath the peaks between 3.2 and 3.4 /im, 8.5 and 10 
^m, 11 and 12 ^m, etc. 



of empirical data optimized on molecules which necessarily differ from those 
of interest here. Thus, a difference of ^ 0.085/im was noted between the 
IR intensity spectrum of Fig. 10 obtained with the semi-empirical method 
PM3 (center graph) and the emission spectrum (bottom graph) obtained with 
the MM code. In the latter figure, the spectrum was shifted in the opposite 
direction to restore agreement with the more accurate PM3 method. 

As an aside, it must be stressed that the broadening demonstrated in Fig. 8 
and 10, as well as the observed UIB widths, cannot be approached by any 
reasonable thermal Doppler broadening, as the velocity of heavy molecules in 
the cold astronomical environments is a tiny fraction of the velocity of light. 
Broadening by rotational bands must also be negligible, because of weight, 
cold and rarity of collisions. 



20 R. Papoular 



4 A TENTATIVE UIB-LIKE SYNTHETIC 
SPECTRUM 



Even allowing for line broadening, the line spectra of intermediate sized molecules 
(100-500 atoms) are not congested enough to yield continuous spectra resem- 
bling observed UIB spectra. One way forward is therefore to increase the 
molecular size. However, as stated above, a serious treatment of molecular 
emission becomes intractable beyond a size of 100 atoms or so. An alternative 
is to increase the number of tractable molecules and add their computed spec- 
tra. In this direction, an informed choice of adequate molecules can be guided 
by a study of the integrated absorption coefficients, /(i^), of a large number 
of CHONS molecules, which was done by Papoular p011b[) . The structures 
that were considered fall in general categories: aliphatic chains, aromatics, 
coronene-like, trios (a pentagon between 2 hexagons) and their associations. 
It was possible to associate each category with a characteristic congregation of 
absorption lines. Building upon this data, 21 structures were selected, ranging 
from 27 to 112 atoms in size, and their emission spectra calculated according 
to the procedure described above. Each was ascribed an adjustable weight, 
which was determined by cut-and-try so as to steer their weighted sum towards 
the shape of an arbitrary but characteristic UIB spectrum: that of NGC1482 
(Smith et al. 2007). A preliminary result is drawn in Fig. 11. 

Even though no line broadening is involved here other than that due to finite 
time FFT (a negligible 0.0035 /xm), and that which is due to anharmonicity 
and computed by the semi-empirical code, as explained above, the spectral 
intensity never falls to zero beyond 5 /im. Also, as expected, several bands 
are formed where lines congregate in the absorption spectra of the molecules. 
It must be stressed again that this result was obtained in the absence of 
any thermal, rotational or pressure broadening, or convolution of the lines 
with any arbitrary Gaussian or Lorentzian profile whose wings would produce 
similar effects. 

Improving the similarity of synthetic to observed UIB spectra requires the 
multiplication of molecules involved, as well as finer tailoring of their spectra 
by introducing small changes in their structures and in the details of their 
excitation. This is postponed to later work, as it falls outside the scope of the 
stated purpose of this paper. 
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Figure 11. In black: the calculated emission spectrum, E{v), of a collection of 22 CHONS structures described in the text 
and totalling 1110 atoms (3324 vibrational modes); the area under the curve is filled. Red line: for comparison, the measured 
flux from NGC1482 (Smith et al. 2007). Magenta dots: the detected features listed by Smith et al. (Table 3), plus the CH 
stretching features at 3.29 and 3.4 /^m. All 3 sets of ordinates are arbitrary. Note the budding continuum between 6 and 14 
/im, and the plateau between 16.5 and 18.5 /xm. No line broadening is involved other than that which is due to anharmonicity 
and computed by the semi-empirical modeling code. 



5 FROM MOLECULES TO GRAINS 



At this stage, it has become clear that the synthesis of the various observed 
spectral emission bands and plateaus requires a very large number of medium- 
sized molecules or a smaller number of larger molecules, such as may provide 
the necessary diversity of functional groups, accompanied by dense enough 
spectral line densities for the lines to start merging together under reciprocal 
attraction due to anharmonicity. On the other hand, if the observed underlying 
continuum is due to structures of the same composition as the said molecules, 
then these structures must be still larger. That large enough molecules exist 
is made possible by collisions in the ISM. Depending on the C/0 ratio of the 
environment, silicate or carbonaceous grains will form preferentially. Both are 
subject to the mechanisms described above for the formation of bands (see 
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Fig. 7 to 10). This evolution can be viewed as taking place in two steps. 
First, upon chance encounters, molecules can stick together through van der 
Waals attraction and H bonds, thus forming weakly coupled clusters. As these 
grow larger with time, they become more resilient and covalent bonds may 
ultimately be formed within, giving finally birth to grains. 
How do the relative contributions of bands and continuum to dust emission 
vary as the grain size increases? Unfortunately, present common modeling 
packages can hardly handle clusters with thousands of atoms. It is possible, 
however, to forsee some of their general properties. Using the capabilities of 
available chemical simulation and visualization codes, one observes graphically 
that normal modes associated with functional groups flock below ^ 20 /im 
while modes above are associated with motions of a large part, or the whole, 
of the structure. The latter are called bulk or skeletal modes and often show 
up as "slow" deformations of the structure, which can be identified with the 
phonons of the solid state (so-called "drumhead" modes). As the molecular 
size increases, phonons of proportionally longer wavelength are created, in- 
vading the far IR spectrum where they form a continuum. 
According to Debye (see Rosenberg 1988), the upper limit of the phonon 
frequency in crystals is 

tJph = (6^2^)l/3i, , 

where A^ is the atomic density of the solid, and f , the sound velocity. Typically, 
on Earth, N = 10^'' at.cm^"^, and v = 2000 m.s^^, leading to a Debye wave- 
length of ^ 20 /im, defining a lower limit to the photon spectrum. Moreover, 
Debye's elementary theory also shows that the spectral density of phonons 
increases as the square of their frequency, i.e. maximum density occurs at 
Debye's wavelength. This is at the origin of the frequent assumption that the 
far IR emissivity of grains goes like A^^ (see, for instance, Zhang et al. 2011). 
It must be kept in mind, however, that the UIB carriers cannot be as densely 
packed as earthly crystals are, for this would impede the formation of func- 
tional groups. So, Debye's limit in IS grains likely lies somewhat beyond 20 
/im, as suggested by the spectral energy distributions emanating from circum- 
stellar shells. But, in spite of a lower density, the relative contribution of grains 
to the near and mid-IR is bound to decrease as their size increases, for the 
phonons claim a growing fraction of the excitation energy: for the 493-atom 
molecule, this fraction is already 30 %. 

At the lower end of the size distribution, the contribution to UIB emission is 
limited by the condition that the structure of the functional groups not be 
perturbed severely enough by one exciting event, as to change the spectrum 
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of their normal mode intensities beyond recognition. Let us arbitrarily define 
this limit by a maximum of two quanta per mode upon excitation. Assuming 
only half of the modes are effectively excited (see discussion in Sec. 2), and 
the weakest quantum energy to be 0.1 eV (A ^ 13/im), the limit is defined by 

where E is the excitation energy in eV. For excitation by H capture, E — 4.5 
eV; for UV excitation, E — 20 eV. The corresponding size "cut-offs" are 17 
and 69 atoms, respectively. Of course, these are to be considered only as orders 
of magnitude. 

There is also the case of grains in thermal equilibrium with a luminous star in 
the vicinity. For such grains to be detected in the mid-IR, their steady-state 
temperature must be about 300 K. In this case, the average energy per atom, 
E{eV)/N — T{K)/7800, is too small to be of concern, whatever the grain 
size. 



6 BANDS IN ABSORPTION 

One may wonder how our medium and large molecules behave in absorption. 
Unless they are maintained somehow in a state of permanent excitation, the 
molecules will remain in the ground state most of the time. It was stated in Sec. 
2.1 that molecules in the ground state carried only discrete lines because the 
anharmonicity mechanisms are not active. Such lines can only be detected in 
absorption if their carriers are hugely concentrated, which is not the case in CS 
shells, and still less in the ISM. Apparently, no such line has yet been observed. 
The core of the difficulty probably lies in the stellar background continuum 
against which the absorption line must be observed. The fluctuations, AB, 
of this background, B, set a lower limit to the detectable absorption signal. 
Indeed, let Siy be the spectral resolution of the detector, and Aj/ the equivalent 
width of the line to be detected. By definition of the latter, 
ABSi' = BAv. 

But Av can also be expressed, in cm~^, by 
Au = j^adv=^, 

where <j is the absorption cross-section in cm^, / the IR intensity in kni/mol, 
and N, the line density of absorbers along the line of sight, in cm^^. From 
this, the minimum line density of absorbers is deduced as 

N^610-f^. (5) 

For instance, if ^ = 0.01 ,5v = 10 cm"-'^ and / = 100 km/mol, then N = 
6 10^"'' cm^^, which is hardly plausible in a CS envelope. 
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For carriers to be observable in absorption, they have to be much larger still 
than those studied above, definitely the size and density of solid grains. For 
such objects, the functional groups that carry the UIBs will be destroyed or 
their features will be drowned in the continuum. This seems to be the case 
at least for the CH aromatic stretching band which, to my knowledge, was 
never observed in absorption; instead, one observes a band extending from 
^ 3.3 to ^ 3.55 /ini, characteristic of aliphatic CH stretchings in methine 
(CH), methylene (CH2) and methyl (CH3), both symmetric and asymmetric 
(see, for instance, Pendleton et al. 1994, Brooke et al. 1999). Figure 10 above 
suggests a tendency towards this direction of change. This is favored by the 
formation of dense internal webs of essentially aliphatic carbon structures. 



7 SUMMARY AND CONCLUSIONS 

In the laboratory, the IR absorption spectrum is the easiest route to under- 
standing the structure of molecules; hence the universal use of IR intensity 
spectra^ I{v)^ each line of which corresponds to one of the normal modes of 
vibration. In space, detecting absorption lines requires unusually large con- 
centrations of material; emission lines are much more common. It was shown 
above that absorption spectra are poor predictors of emission spectra. Not 
unexpectedly, they must be supplemented by the knowledge of the energy 
content of each vibration mode of the emitter. 

If the various normal modes remained isolated and strictly harmonic even in 
the excited state, then, each one would emit the whole energy that was de- 
posited in it in the first place, regardless of its IR intensity. However, as soon 
as the molecule is energized and no longer in its ground state (in which it 
is, in the laboratory), anharmonicity sets in, and induces internal couplings 
between modes, which help redistribute the deposited energy. The equilibrium 
energy distribution among modes depends on the particular molecular excita- 
tion process (thermal, single UV photon absorption, capture of an H radical, 
etc.), but mostly on the internal couplings. Thus, for each candidate carrier 
of an astronomical spectrum, one cannot dispense with the computation of 
the dynamical equilibrium reached shortly after the deposition of excitation 
energy. The analysis is made easier by the fact that the intramolecular energy 
redistribution is much faster than the IR emission process. 
During the emission process, the most IR active vibrations are energy-depleted 
faster than the others, which, in turn, refuel the former faster than they them- 
selves radiate their own energy. In the end, it turns out (see Appendix), that 
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a better approximation to the band emission scales with the product of initial 
energy content times the line intensity, times the frequency squared. 

Another effect of couplings between modes is to induce line broadenings and 
shifts. This is incompatible with the use of IR line intensities, which are S 
functions at the normal mode frequencies. A way around this obstacle is the 
recourse to the computation of the time variations of the total electric dipole 
moment of the molecule. The Fourier analysis of these variations deliver the 
emission spectrum. The Appendix shows that the emission at frequency v 
scales with the product of initial energy content times v^ , times the square of 
the dipole moment component at frequency i> (eq. 4). This is proposed here 
as the more accurate of the three estimates; it comes, of course, at a price: 
more elaborate and time consuming computational codes. 

For clarity, the three approaches are first illustrated, here, by chemical com- 
putations on a small hydrocarbon molecule. This is barely enough to unveil 
some essential differences between them, as well as a small line broadening. 
However, our main purpose is to study the formation of bands and plateaus, 
starting from the IR lines. This requires larger and larger molecules. There- 
fore, the dipole moment procedure was then applied to molecules containing, 
respectively, 101 and 493 atoms of the CHON family. Because of the higher 
density of IR lines, many broadened lines merge together, giving rise to fewer, 
but wider, spectral features which are precursors of bands. Moreover, the 
wings of neighboring lines overlap and form plateaus, where line densities are 
highest. No thermal broadening is involved. 

Bands and plateaus can also be obtained by associating a large number of 
smaller molecules. A very rich diversity of intermediate-size molecules can 
produce enough intercalated broadened lines, packed within a few character- 
istic spectral regions, for their sum to look like a continuous band spectrum. 
This procedure is illustrated here by a spectrum resulting from the summa- 
tion of the spectra, computed according to the dipole moment method, of 22 
different CHONS structures, ranging from 27 to 112 atoms in size and belong- 
ing to 5 families of structures (PAH, aromatics, aliphatic chains, etc.). The 
time required for the computation of one structure on a desk-top PC is a few 
hours. 

Beyond about 20 /xm, large molecules exhibit far IR bands associated with 
bulk vibration modes (phonons) involving motions of the whole, or large parts, 
of the molecular structure rather than only the few atoms comprised in the 
functional groups which carry the mid-IR bands. As the molecule grows larger. 
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it emits more and more in these bands and, relatively, less and less in the UIB 
bands. 

These general computational results could not be obtained on the basis of 
the usual normal mode analysis, although the latter delivers essential guiding 
informations in the quest for candidate models of the UIB carriers. 
The so-called semi- empirical dynamics code used here takes care of the exci- 
tation process and of the anharmonicity of the motion, which depends upon 
the molecule, the excitation energy and the particular vibration mode con- 
sidered. If the excitation energy is high enough, overtone and combination 
bands automatically emerge from the analysis of the variations of the dipole 
moment in time. This analysis delivers the emission without going through 
the Normal Modes. No assumptions are added to shift or broaden the bands 
artificially. Plateaus arise naturally from overlapping band wings. Each of the 
strong UIBs beyond 5 /xm can be associated with one family of structures. 
While DFT (Density Functional Theory) and its derivatives also deliver accu- 
rate Normal Mode intensities and frequencies, they are basically static meth- 
ods, so they are not directly applicable to analysis of dipole moment variations. 
To my knowledge, even Time-Dependent DFT (TDDFT) has been applied 
only to adiabatic transitions between equilibrium states, not to dynam,ics. 



8 APPENDIX: RADIATED ENERGY 

Present chemistry simulation codes do not include spontaneous radiation. A 
simple analytical argument can, however, be made to deliver an operational 
approximation to reality. 

When energy is deposited in a molecule, several modes are excited; which ones 
depends on the excitation process. If the modes were strictly harmonic and, 
therefore, completely uncoupled, each and everyone would start radiating in 
the IR, at its own rate, different than the others. Assuming the excitation rate 
to be slow enough, each mode would ultimately radiate the energy which was 
initially deposited in it; the radiated power would be obtained by multiplying 
this with the frequency of excitation events. But, in a real molecule above its 
ground state, anharmonicity cannot be overlooked and the modes are coupled 
more or less. As a consequence, intramolecular vibration redistribution then 
quickly changes the initial distribution into a new one, which is mostly in- 
dependent of the excitation process, but characteristic of the molecule itself. 
For simplicity, assume first that the molecule has only 2 vibration modes, and 
let this characteristic distribution run like cj' , with i = 1,2, where q is a 
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pure, constant number. Let Ei be the energy, t^ the radiative lifetime of the 
ith mode, and e the energy exchanged between modes. Then, 

Ei:^-Ei/Ti+e (6) 

E2 = -E2/T2 - e (7) 

e = a{ciEi - C2E2), (8) 

where the dot designates the derivative and a is the mode couphng constant. 
Note that, when Ei oc c^ , energy exchange stops. It is readily verified that 
the solution of this set is of the form 

El — aiexp{—t/T) + (3iexp{—t/a) (9) 

E2 = a2exp{-t/T) + (32exp{-t/a) (10) 

e = aexp{—t/T) + l3exp{—t/a). (11) 

Identifying separately the coefficients of exponentials, and applying initial 
conditions, it is found that 

(12) 
(13) 



(14) 
(15) 
(16) 
(17) 

" aici + C2) ' ('^^ 

where the subscript Oi designates the initial value of the mode energy. 
The meaning of these equations is clear when, as is usually the case, the 
coupling constant, a, is strong enough that cr <^ t. In a first, short, phase (the 
transient, with time constant cr), radiation is insignificant, but the energies in 
the two modes are forced to quickly readjust so as to reduce e (eq. 8). They 
then can relax at the same rate, r, intermediate between the two radiation 
rates ti and T2, while energy is being continuously transfered from the less 
radiative to the more radiative mode, at a very slow rate. 
Finally, the energies radiated by modes 1 and 2 are, respectively 

Eidt/ri ~ tui/ti = T2Q!i (19) 

ClTl + C2T2 

/ E2dt/T2 = Ta2/T2 = — Tia2 , (20) 
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whose sum equals ai + a2 — Eqi +Eq2, as it should, since there is no outlet for 
the initial energies other than radiation. The ratio of total energies radiated 
at the respective mode frequencies is 

^^. (21) 

a2 Ti 

The forms of the solutions above are easily generalized to any number of 
modes, assuming a uniform coupling constant. The coupling constant only 
impacts the transient time constant, a, which remains very short anyway (a 
few tens of picoseconds) . One can therefore identify Eoi with aoi , the charac- 
teristic energy distribution; the ratio of radiated energies then become 

fl^. (22) 

It is useful to relate the radiative time constants to the IR intensities delivered 
by the chemistry codes, so as to expedite the computation of the emission spec- 
trum. Indeed, the radiative lifetime is the inverse of Einstein's spontaneous 
emission coefficient (see Wilson et al. 1955) 

A{y) = ^)J^"^(^) ' (23) 

in electrostatic CGS units, and written in simplified form, for one particular 
transition and overlooking degeneracies. Here, miy) is the dipole moment 
oscillating component at frequency v. Now, it comes as no surprise (see Atkins 
1997) that the corresponding IR intensity is also linked to m(v) by 

I{v) = 2.%vm{vf , (24) 

where m is in Debye and v in cm~^. Thus, Einstein's coefficient scales like 
v^Iiy) and v^mii/f' . Finally, the energy ultimately radiated by mode v scales 
like 

E* <^v'^m{vfEQ(y), 

where E^iy) is the energy initially deposited in that mode (that is, immedi- 
ately after the very quick redistribution of the total energy provided by the 
exciting process). This is the expression used in Sec. 2.3. It is a correction to 
the presentation of the dipole moment spectrum proposed in Papoular pOOip . 
It should be noted that Einstein's coefficients are generally unknown for the 
proposed candidate carriers, and their determination is far from trivial. Pre- 
vious authors in this field have turned around the difficulty by pointing out 
their linear relation with the corresponding absorption cross-sections. Even 
those are very scarce, having been obtained by measurements on a few partic- 
ular molecules. By contrast, the procedure outlined here is easily applicable 
to any structure, provided it can be handled by the computer at hand. 
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